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The production of antihydrogen in flight in p-nucleus collisions is calculated theoretically in the 
Plane Wave Born Approximation (which is equivalent to the straight line semiclassical approxima- 
tion). Antihydrogen has been produced in this way at LEAR/CERN ||^ and is presently studied at 
FERMILAB at various p energies. Dirac wave functions for the leptons are used, taking first order 
{Za) corrections into account. Analytical results are obtained for differential cross-sections. Total 



^\ ' cross sections are obtained by numerical integration. The dependence on the transverse momentum 



transfer is studied and the accuracy of the equivalent photon approximation and a recent variant 
by Munger, Brodsky,and Schmidt [Q| is discussed as a function of beam energy. 

PACS numbers: 25.43, 32.10, 34.80, 41.20 



Csl ■ I. INTRODUCTION 

> : 

. Antihydrogen in flight was first produced and detected at CERN/LEAR using the process 

(N ; 

p + Z^H + e^+Z, (1.1) 

Qs^ ■ where Z=54 (Xe). This process was first suggested by Munger et al. ||], where also calculations of the cross section 
. for the antihydrogen production were performed in the equivalent photon approximation (EPA). A similar pair 
production process with e-capture will occur with large cross section at relativistic heavy ion colliders Q. This is 
Oh important because it leads to a beam loss. Theoretical calculations were done by many groups |^-^. They were 
mainly interested in the limit of high energies. The applicability of EPA will become more and more questionable for 
^ ' the lower energies, where the H production experiments were done |l]j^. The dominant graph is shown in figure 1. It 

' was shown in that bremsstrahlung pair production with capture can be safely neglected. 
j_j ■ By crossing symmetry the process 7 + p — > e~ -I- H is related to the photoeffect 7-I- H — > P + e~ [§) which is 
1 treated in the literature (see e.g. j|]). These calculations provide approximate theoretical values for the H production 
cross section. This is also relevant for the ongoing experiments at Fermilab , where also the energy dependence of 
the process will be investigated and possible future experiments with fast H beams, like Lamb shift measurements 
On the other hand, we have to investigate the accuracy and limitations of the EPA. It is the purpose of this 
paper to carry out an accurate calculation of the cross section for eq. (1.1). Analytical expressions are given. 

II. PRODUCTION OF H WITH ANTIPROTON BEAMS 

In the semiclassical approach we assume that the relativistic antiproton moves along a straight line. In its frame of 
reference the time-dependent electromagnetic fields of the target nucleus is given by 



*Dedicated to the memory of Franz Osterfeld. 
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A(r, i) = V 0(r, t), where (p{r,t) 



Zre 



(2.1) 



\r-r'{t)\ ' 

where r = {x,y,z), and r' — {bx,by,"/vt), where v is the relative velocity, 7 = (1 — v^)^^/^ , and Ztb is the target 
charge. We can write cj) in the integral representation 

giq.[r-r'(t)] 



(r,t) 



27r2 J - q 

In the first-order perturbation theory, the H-production amphtude for a colhsion with impact parameter b 
+ 6^ is given in terms of the transition density and current, p (r) and j(r), respectively, by 



(2.2) 



ai«t = i j dte'^' j d\ [p(r)0(r,O-j(r).A(r,i)] 

giq.[r-r'(t)] 



dt e" 



d-^q / d r 



[P (r) - v.j(r) 



(2.3) 
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where uj — e + m is the sum of the electron, e, and the positron energy, m (the binding energy of H can be neglected). 
In the last step we have used eqs. (2) and (3). Performing the integral over time and using the continuity equation, 
V.j = —iiop, we get 



Ztg 



d^r e'^ " 



jz{r) , qt.j 



+ 



(2.4) 
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where the index z (t) denotes the components along (transverse to) the beam direction, and Q = (qt,a;/v). The 
transition current is given in terms of the Dirac matrices, 7?, and the lepton wavefunctions j(r) = e [5'*^^-'] l^^^^K 
Thus we can rewrite the H-production as 

Zre r 6*1* '^ 



d'qt 



F{Ct), 



(2.5) 



where 



(r)e^Q- (^-H^i^)vI/(+)(r) 



ie K . / d^r 



(r) e'^ '' 7° 7^ *(+)(r) 



(2.6) 



In the last equality we used the definition K = (Kj, K^) = (qt/ij-', l/v7^) and (7°, ^) arc the Dirac matrices, where 
we follow the procedure of rcf. ||^. The total cross section is obtained by integrating the square of the expression (8) 
over all possible impact parameters: 
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d^t 
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(Q2-^2)S 



(2.7) 



spins ^ ' spins 

The same result can be obtained in the Plane- Wave Born approximation, as shown in the appendix C. 

Using the positron and the electron wavefunctions as given by eqs. (A. 4) of appendix A and cq. (B.l) of the 
appendix B, we get 



F(Q) ^iej d^r {u e-'P " + ( t^.K) e''^-'- 

2m J 
ie I d^r Jue'('5-P) '-(7^.K) 



1 + ^/^.^) 
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where in the last equation we neglected terms of highest order in Za. 
Integrating by parts 



i^(Q) = ie y rfV {u e'(Q-P)- ( T^.K) (l + ^ 7°^-^j v*„„„_,(r) 
+ ■Fe*C!-(l^.K)vvE'„o„_,(r)} 



l + i/7^-(Q-p))v(e 



-r/oo 



Q-P 



(2.9) 



Since the corrections are in first order in Za, we have replaced ^non-r{s) by its constant value Xj^aJ"^ in the last 
term of eq. (2.9). If we would do the same in the first term it would be identically zero (except for Q = p). This is 
the reason why we need to keep the corrections in the wavefunctions to first order in Za\ for e ^ m these corrections 
yield a term of the same order as the plane- wave to the total cross section, as we shall see later. We rewrite eq. (2.9) 

as 



5/2 u^v 



where 



with 



F(Q) = 4ieV^ao , 

(Q-p) 

A = a(l^.K) + (1^.K)7° {-^.h) + (l^.d)7°(l^.K) 
1 e ^ Q-p 



(2.10) 



(2.11) 



c = — - 



(Q-p)' Q^-P^' 2m (Q-p)' 

Q-P 
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2m(Q2_p2) • 

Inserting (2.10) into (2.7), the cross section becomes 

_ 16e2 /ZtcV f (fp f 
""'^0 \~) J (2^io ^^2_„,2^2 
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spins 



The sum runs over the spins of the positron and the electron. Using p^dp = spds, and standard trace techniques 
we get 



da 32e2p f Zre 



d% 



1 



iQ'-co^r (Q-p)' 
X I (e - 1) (b - d)' + 4 (K.b) (K.d) 

+ (e + 1) a^K^ + 2a [2 (p.K) (b.K) - p.(b - d)K'^] 



(2.14) 



Using the definition of K, it is illustrative to separate the integrand of the above equation in terms of longitudinal 
and transverse components: 

2 " 
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where 
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(2.17) 



In these equations e = qt/gt is a unit vector in the transverse direction. The cross sections cr-yT, C7*l, and a^*LT 
are interpreted as the cross sections for H by virtual transverse and longitudinal photons and a interference term, 
respectively. Note that the transverse and longitudinal directions are with respect to the beam aa;zs,nor with respect 
to the photon momentum, as is usually meant by this term. Only for 7^1, this definition agrees with the definition 
of transverse and longitudinal virtual photons. However, this separation is very useful, as we will see next. 



III. PRODUCTION OF H BY REAL PHOTONS 



The cross section for production of H by real photons is given by [g[ 

da^ = 2Ti\Vf.,\^ 5{uj-e-l)-'^^^ 



{2^r 



where 



— / d r 

UJ 



^(-)(i.)l e*'"-'" (e.-^)^'W(r) , 



(3.1) 



(3.2) 



with K = 2,Lo/c, where z is the unit vector along the photon incident direction, and e is photon polarization unit 
vector. 

Using the positron and the electron wavefunctions as given by eqs. (A. 4) of appendix A and eq. (B.l) of the 
appendix B, and performing similar steps as in section 2, we get 



87rV2 uBy 
,5/2^1/2 (,_p)2 ' 



(3.3) 



where 



B = a,(l^.T?) + ( 7^. K*)7° ( T^.b.) + ( 7^.d„)7"( 7^.^) 



(3.4) 



where a^, b^, and are the quantities defined in eqs. (2.12), but with Q replaced by 7? —'iuj/c. Inserting these 
results in eq. (3.1) and summing over spins we get 
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(k-p) 



(e-l) 



(b,-d«y - 



(p.S)^ 



(e + 1) + 2a«, 



(3.5) 



We notice that the above equation can also be obtained from eqs. (2.15-2.18) in the hmit — > and v — > c. In 
this hmit a-y*L — > 0, and ay-LT — >■ 0, and aj-T becomes the cross section for the production H by real photons. 
Integrating eq. (3.5) over the azimuthal angle we get (without any further approximations!) 



27r 
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{s — p cos t 



e — p cos 9 



Lu{e + 2) 



(3.6) 



In figure 2 we plot the cross section for the angular distribution of the electrons in H production by real photons. 
We observe that the higher the electron energy is, the more forward peaked the distribution becomes. As we can 
immediately deduce from eq. (3.6), for e ^ 1, the width of the peak is given by A cos 6 ~ 1. For e ~ 1 the distribution 
is proportional to sin^ 9. 

Integrating eq. (3.6) over 9 we get 

47re^ 



cr (7 + p 



H 



P 



^ +r+3 



g + 2 
P 



In (e + p) 
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It is instructive to derive this cross section by a different method. In some textbooks (see, e.g., ref. 0,^) one can 
find a calculation of the cross sections for the annihilation of a positron with an electron in the K-shell of a nucleus. 
From charge conjugation invariance (which is valid for QED), this cross section is the same as the cross section for 
the annihilation of an electron with a positron in the K-shell of an H. The cross section is |To| , pl| 



H 



27re2 1 



7 + P) 



2 4 

r+3 



p 



■ In (e -1- p) 
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Using the detailed balance theorem, the above equation yields exactly the same cross section of H production by 
real photons as given by eq. (3.7). The reason is that the calculation leading to eq. (3.7), originally due to Sauter 
p^ , used the same corrections to Za order for the electron and positron wavefunctions, as we did to obtain eq. (3.1) 
(see Appendices A and B). 

The relevance of these corrections can be better understood by using the high energy limit, e 3> 1. From eq. (3.7) 
we get 

47re2 1 
an £ 



o-(7 + P 



H 



(3.9) 



If in our calculation, we had used plane-waves for the electron and the hydrogenic wave function for the positron, 
without the corrections to order Za, the cross section for the production of H by real photons would be given by 

da^ 

■(LJ) =4:n-^- 

PWA "o^ [e-pcosi 



d9 ^ ' 



•,4 ' 



(3.10) 



the integral of which being 



"fiPWA 



(3.11) 



In the high energy limit: 
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\^i\pwA ~ ^ (3-12) 

We see that the difference between this resuh and that of eq. (3.9) is a factor 8/3. This is just the reason why 
the corrections of order Za have to be inciuded in the caicuiation of H. They yield terms to the matrix elements for 
photo-production of H of the same order as the terms of lowest order. The origin of the difference between the two 
results are the small distances which enter in the calculation of the matrix elements of eq. (2.9). The corrections 
are essential to account for their effects properly. We note that these corrections are enough to account for a good 
description of H production. 

The connection between the production of H by real photons and by virtual ones will prove to be very useful, as 
we shall see in the next section. 



IV. PRODUCTION OF H BY VIRTUAL PHOTONS 

As we have seen in last section, the production of H in a collision of p with nuclear targets is directly related to 
the production cross section by real photons in the limit — > and v — > c. It is important to determine at which 
values of qt and 7 this condition applies. In figures 3 and 4 we plot the production cross sections by longitudinal 
and transverse virtual photons, cr-yT, cr-^.^, respectively. These are the integrals of eqs. (2.16) and (2.17), which are 
performed numerically. The integral of the interference term, cr-y'LT: eq. (2.18), is very small and is not shown. In 
figure 3 the cross sections are calculated for e = 2, and 7 = 3, and in figure 4 for 7 = 10 and e — 2. We observe that 
for low 7's the cross section cr-y^L, is as important as u-^-t-, while for large 7's the cross section ct^.t dominates. For 
large 7's the cross section a^-j^ is only relevant at qt ~ 0. But this region contributes little to the total cross section, 
which is given by 

^7 = — — de d{ql) ^* [(j^>T{Q:U:) + cj^.L{Q:U}) + a^.L{Q,u;)] , (4.1) 

TTUJ \ y J Ji Jo {Q^ - uj^Y 

since the first term inside the integral suppresses ~ values. This can also be seen by calculating the total cross 
section for H production as a function of 7. This is shown in figure 5 where the solid curve shows the cross section as 
given by eq. (4.1), and the dashed line is the cross section with only the inclusion of a-y-T- We see that for 7 of the 
order of 8, or larger, the cross section is dominated by the transverse virtual photon component. Moreover we also 
observe that the cross section for large 7's disagree with the equivalent photon approximation (EPA) used in ref. |^ 
to calculate the cross section for the production of H. It is thus convenient to study under what circumstances the 
equivalent photon approximation is reliable. 

It might appear strange that the longitudinal part of the cross section is relevant until such large values of 7. It is 
thus reasonable to check this result with a more schematic calculation. This can be achieved by using plane-waves 
for the electron and hydrogenic waves for the positron, without the correction terms to order Za. In this case, we get 



^(Q) . 4.e0Fa-/^ , 

(Q-p) 



which yields the cross section 



da _ 1 (ZreV [°° ^ qf 



(4.2) 



(4.3) 



dedn TTUj \ Y J Jo ' (Q2 - u'^f 
where the interference term is exactly zero, and 

da^ 32e^p 1 ^da^ 32fpu^ f 1 V 

— (Q,ll!) = — r o , and — (Q,^^) = ^-^5— ^5- s ■ (4-4) 
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The integral over in the expressions above can be done analytically yielding 
(t^.t(Q,w) 1 _ 128^6^^ (p' + 3Q2) (3p2 + fl 

Using this result, the integral over in (4.3) can also be performed analytically but resulting in too long expressions 
to be transcribed here. In figure 6 we plot 1 — {da l / de) / {dax / de) obtained in this approximation, for e = 1, 2, and 
3, respectively, and as a function of 7. We observe indeed that even for relatively large values of 7 (e.g., 7^6) the 
longitudinal part of the cross section is still substantially relevant for the calculation of the total cross section. 




V. THE EQUIVALENT PHOTON APPROXIMATION 

The equivalent photon approximation is a well known method to obtain cross sections for virtual photon processes 
in QED. It is described in several textbooks (see, e.g., It is valid for large 7's, in which case the cross section is 
dominated by transverse photons, as with the production cross section of H. In figure 7 we plot the virtual photon 
cross section a-y^T {Q, ^) for £ = 1,5, and 10, for 7 = 100, and as a function of qt- In this logarithmic plot it is evident 
that the function a~f*T {Qi^) is approximately constant until a cutoff value, q™^^, at which it drops rapidly to zero. 
On the other hand, the function 

which we call by "equivalent photon number", is strongly peaked at qt ~ W/7V, which is very small for 7^1. Thus 
it is fair to write 

d(j 1 1 

— d{ql) -n{^,Q) a-,,T{qt=0,u:)^-Oy{uj) d{ql)n{^,Q) (5.2) 

as Jq oj Lo Jq 

where in the last equality we have approximated ct^t (qt = 0,a;) ~ (w) . This approximation is indeed valid for 
7 ^ 1, as we can see from figure 8, where we plot the ratio between these two quantities for e = 1, 5, and 10, and 
as a function of 7. For 7 bigger than 10 the approximation is quite good, and it is even better for the lower values 
of e. Eq. (5.2) is known as the "equivalent photon approximation". The problem with the approximation is that the 
integral in eq. (5.2) diverges logarithmically. The approximation is only valid if we include a cutoff parameter g™'*'^, 
determined by the value of qt at which a-yT {Q, ^) drops to zero. A hint to obtain this parameter is to look at figure 
7. We see that the transverse virtual photon cross section drops to zero at qt ~ e. Another hint comes from figure 
9, where we plot the energy spectrum of the electron, da/ds, for 7 = 3, 10, and 10, respectively, obtained from a 
numerical integration of eq. (4.1). We observe that the energy spectrum peaks at e ~ 2 — 3, irrespective of the value 
of 7. We thus conclude that an appropriate value of the cutoff parameter is g™^'^ = 2. Inserting this value in eq. (5.2) 
we get 



da 1 (ZreV 1 , . 

— = - - CT^ (w) 

de TT \ V / LO 



In {Ax + 1) - 4 



4a; + 1 



„2 



where x — — — 5 — . (5-3) 



The dotted curve shown in figure 5 is the integral of the above equation over lu. We see that this approximation 
is very reasonable, for large 7. It is only for the very extremely large 7's that we obtain cr = (2.86 pb) Zi^ In 7, 
where the number inside parenthesis comes from {2a/T:) dm {lu) /uj — 2.86 pb. Only at 7 10'^, and larger, 
the approximation presented in ref. [|j becomes close (in the range of 10%) to the result obtained with the equation 
above. 
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VI. CONCLUSIONS 



We have calculated the cross sections for the production of antihydrogen in collisions of antiprotons with heavy 
nuclear targets The cross section is well described either by a semiclassical method or by the plane-wave Born 
approximation. In fact, both approaches yield the same result for relativistic antiprotons. The cross section can be 
separated into longitudinal and transverse components, corresponding to the velocity of the incident particle. For high 
energies, 7^1, this corresponds to the usual meaning longitudinal and transverse components of the virtual photon. 
A very transparent and simple formulation is obtained using the lowest order corrections to the positron and electron 
wavefunctions. At ultra-relativistic antiproton energies, the contribution of the longitudinal virtual photons to the 
total cross section vanishes. The remaining terms of the cross section can be factorized in terms of a virtual photon 
spectra and the cross section induced by real photons. However, the factorization depends on a cutoff parameter, 
which is not well defined. This leads to differences up to a factor of 2 between our results and those based on the 
equivalent photon approximation. Production of antiatoms with the electron at low energies, typically e ~ 2, is 
favored. The angular distribution of the electrons is forward peaked within a angular interval ~ 1/e. 

The equivalent photon approximation is only of help as a qualitative guidance for low values of 7 (7 smaller than 
about 10). This is the case for the CERN/LEAR experiment {p = 1.94GeV/c), to a lesser extent for the Fermilab 
experiment Q where the beam momentum can vary from 3.45 to 8.8 GeV/c. The formulation presented here allows 
for a better quantitative calculation of the production of antiatoms at these energies. The discrepancy to earlier 
calculations for small values of 7 is at present not understood. 
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APPENDIX A: POSITRON WAVEFUNCTION 

Here we deduce a first order (Za) correction of the positron wavefunction, important in the calculation of the H 
production cross section. To be general, we consider the production of a bound positron in a relativistic antinucleus, 
with charge —Z, incident on a target nucleus with charge Z^- Much of our calculation is based on arguments presented 
in ref. |^ with connection to the photoelectric effect. 

The Dirac equation for a positron in the field of the antinucleus is: 

ex]/ = ["^.(p + eA) - /3r7i - e0] . (Al) 

For A = : 



e — U + /3m + 



* = , (A2) 



where U = —e(j> = —Ze^/r is the Coulomb field. The positron will be most likely be produced at the s-state of the 
antinucleus. To lowest order the wavefunction is given by the non-relativistic hydrogenic wavefunction 

"Unon-Ar) - ^ (-) ' 6-'-/'^° , (A3) 

where ao = Xj [Ze^ra). 

To first-order, a corrected wavefunction (to order Ze^) is given by 

where v denotes the positron spinor. 

Applying the operator e — [/ — m/3 + to (A. 2) we get 

(V^ + - 2eU) ^' = {i^^U - U'^^ (A5) 

Substituting (A. 4) into (A. 5) and expanding: 

(V^ + p2 - 2eV) *„o„-rV + (V^ + p2 _ 2e[/) 



(A6) 



Since = + rr? ~ — 2m|es|, "^^^^ oc Ze^, and U oc Ze^, we get to lowest order 
(V' +p2 _ 2eU) *„o„-r ^ (V^ - 2m|e,| - 2|e,|f/) 

and eq. (A. 6) becomes 

(V^ - 2m|e,| - 2toC/)*(i) = a^.^C/) ^'„on-rV , 
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- l^^l - — ^ =-^-[ "^non-rV , (A7) 



or 

The non-relativistic wave function obeys the equation (V^/2m — \ss\ + Ze'^/r)^non-r = 0) from which we deduce 
that 

- \£s\ - —y^^non-r = - —) *non-r , (A8) 

Thus, if ^^^^ = i/{2m)vli ^'^non-ri it will be the solution of (A. 7). An approximate solution of (A.l) is therefore 



^(+) 



1 + — .V 
2m 



1 + ^/^.^ 
2m 



Y^non-r , (A9) 



The relevant distances for the non-relativistic wavefunction are r ~ l/{mZe'^). The correction term should be 
good within these distances. But, for the ground state (or any s-state) it can be used for any value of r, since the 
derivative of the exponential function (A. 3) is always proportional to Ze^. Because of that, we can use the corrected 
wavefunction in our calculation of H production where, as we see in section 2, the small values of r are essential in 
the computation of the matrix elements. 

APPENDIX B: ELECTRON WAVEFUNCTION 

For the electron wavefunction wc use a plane wave and a correction term to account for the distortion due to the 
antinuclcus charge. As in Appendix A, the correction term is considered to be proportional to Ze^. The wavefunction 
is then given by 

= u e'P "" + . (Bl) 

In section 2 we show that only the Fourier transform of ^' will enter the calculation. This Fourier transform can 
be deduced directly from the Dirac equation for the electron in the presence of a Coulomb field of an antinucleus: 



(7°e + i *' = 7°u e^P■^ (B2) 

Applying on both sides of this equation the operator ^7°e + i'^.^+m^ we get 

(A + p2) = Ze^ (7°£ + i T^.^+mj (7°u) ^ . (B3) 
Multiplying by e""! "" and integrating over d^r we get 

Air 

(p2 - q') = [27°£ + i 7^.(q - p)] (7°u)^^--^ , (B4) 
where we have used the identity (j°e + i'^.^—m^ (7°u) = 0. Thus, 



In section 2 we use this equation to calculate the matrix element for the production of anti-atom. 
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APPENDIX C: PLANE- WAVE BORN APPROXIMATION 



In the plane wave Born approximation the transition matrix element is given by 

Tfi = I [pfi (r) .^(r) - j/,(r).A(r)] (CI) 



where 

and (k|r) = e**^ "" is a plane wave for the antiproton. Using 

iai|r-r'| ■. iK.(r-r') 



and Jfi = vpfi, we get 

d'r ^^-^ [pfi (r) - j/,(r).v] = 4nZe (C4) 

where Q = k— k/ and -F(Q) is given by eq. (2.6). The cross section is given by 

da ( eV 



For relativistic antiproton energies 



Ql = ki — kf cosd ki — kf c:^ w/v 

Qt = qt = kf sin {e)c:iEv sin {6) dCl = d'^qt/ {Evf , (C5) 



so that 



spins I Q'^ — up- \ 



The qt integration ranges from to a max;imum value En, where E and v are the antiproton energy and velocity, 
respectively. This value is however much larger than the relevant energies entering the matrix elements in -F(Q). 
Thus, the expression above is the same as the one derived in section 2, eq. (2.7). 

APPENDIX D: FIGURE CAPTIONS 

1. Diagram for H in collisions of p with a nuclear target. 

2. Angular distribution of the electrons in H production by real photons for several electron energies. 

3. Production cross sections by longitudinal and transverse virtual photons, u-^*t, ct^'l, respectively. The cross 
sections are calculated for £ = 2, and 7 = 3. 

4. Same as in figure 3, but for 7 = 10 and £ = 2. 
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5. H production in a collision of an antiproton with a proton target, as a function of 7. The solid curve shows the 
cross section as given by eq. (4.1), and the dashed line is the cross section with only the inclusion of a^*T- The 
dashed line is the result of the approximation given by eq. (5.3). 

6. H production in a collision of an antiproton with a proton target calculated rising a plane-wavefunction for the 
electron and a hydrogenic wavefunction for H. The ratio (daT/de — da l / de) / daT / de is shown as a function of 
7 and for e = 1, 2, and 3, respectively. Eqs. (4.3 - 4.5) were used in the calculation. 

7. Virtual photon cross section cr-y.T {Q-.'jJ) for £ = 1, 5, and 10, for 7 = 100, and as a function of qt- 

8. Ratio between c-yT {qt = 0, w) and a-y (w) for e = 1, 5, and 10, and as a function of 7. 

9. Energy spectrum of the electron, da/de, for 7 = 3, 10, and 10, respectively, obtained from a numerical 
integration of eq. (4.1). 
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